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Abstract 



We calculate the relaxation time self-consistently to study the damping of col- 
lective color modes and the color conductivity in a QGP by deriving self-consistent 
' equations for the damping rates of gluons and quarks to leading order QCD by TFD 

including a chemical potential for quarks. We show that the damping rates are not 
sensitive to the chemical potential whereas color conductivity is enhanced consider- 
ably. 



^ ■ 1 Introduction 

CD : 

J3 ■ Heavy ion experiments at RHIC and SPS have now reached a stage where we can expect 

to probe QCD matter beyond its hadronic state in the quark gluon plasma (QGP) phase. 
For this phase of essentially free colored particles we expect a behaviour similar to the one 
in a QED plasma, that is Debye screening of color electric charges and their interactions, 
as well as collective color modes [|[], j§. 

Collective behaviour might be observable in ultrarelativistic heavy ion collisions as 
long as the lifetime of the QGP exceeds the time scale for dissipative processes. Therefore 
it is important to calculate the damping rate for collective phenomena. 

So far, mostly momentum relaxation rates were considered §§,[§], which allow then 
estimates for viscosity, thermal conductivity and energy transport in the plasma. Some 
work has been done also on color relaxation || , |§ , [0] to calculate the color conductivity 
and the color diffusion coefficients. 

We want to extend that work and calculate the static color conductivity a including 
a chemical potential for quarks. The chemical potential has been estimated to become 
zero only when nuclear transparency sets in and reach values up to 1 - 2 times the 
temperature in the region where stopping occurs || . 

The (infrared) divergence occuring in the calculation of color conductivity, which can 
be related to long range (color) magnetic interactions, will be treated in our approach 
self-consistently by including a finite width v for quarks and gluons which acts as a 
self-consistent damping term. In previous works , || , Q this divergence was regularized 
using a (non-perturbative) magnetic screening mass m g as an additional cutoff parameter. 
We will show explicitely, that a finite width will regularize the divergence self-consistently. 
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The paper will be organized as follows: In section ^ we will get the self-consistent 
equations for the width v of both gluons and quarks to first order perturbation theory 
in the hard thermal loop approximation using the real time formalism of thermal field 
dynamics (TDF) including a chemical potential ^ for quarks. In section || we will discuss 
the self-consistent (numerical) solutions of v and a as functions of temperature and 
chemical potential. We will summarize our results in section 0. 



2 Damping rates of fast moving particles in a QGP 

If one considers a fast moving particle with momentum scale p gT going through a 
weakly coupled QGP (g <C 1), the relativistic dispersion relation for the particle, the 
transverse gluon for instance, can be written according to the poles of its propagator in 
the following form 

E(k) = uj(k) - iu T , (1) 

where i>t is the damping rate of the transverse gluon. 

It is sufficient to consider only the fast moving particles since for the static color 
conductivity the momentum exchange is small (0(<?T)). This results in a damping rate for 
the slow moving particles of order g 2 T in the weak coupling limit, which is small compared 
to the one for the fast moving particles of order g 2 T/ ln(l/g) This corresponds to the 
scenario of a realistic hot QGP, where the momentum of a particle should be of order 
T, that is hard ({p) = 3T), but where the momentum exchange can be rather small. In 
addition, the color diffusion coefficient d c = p ■ u is proportional to the momentum of 
the particle (u^ being the four velocity in the ensemble rest frame), so that color diffusion 
is suppressed for slow particles. 

In covariant gauge the gluon propagator is 

G^ik) = A^{k)A T {k) + B^(k)A L {k) + D^(k)A c {k). (2) 
At and Al are the transverse and longitudinal gluon propagators respectively, given by 

At > l = (3) 

£ is the gauge parameter with associated propagator 

and A^,!}^, D^ u are functions of the loop momentum k, k = (ko, K), 

r fi\ fx k v k o\ k 2 ( k u k 

D,Ak) = (5) 

II is the polarization tensor, and the dispersion relation for the transverse gluon can be 
cast into the form || 




E 2 -P 2 = ^[U]-^l ^F. (6) 
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From Eqs.(jT]) and (||) one obtains the damping rate for a transverse gluon on the mass 
shell g: 

= - 2jF[ ■ (?) 

The crucial step in evaluating Eq.(|^) will be to calculate the imaginary part of the 
polarization tensor. The one-loop graphs for II in SU(N) — QCD are given in Fig]l[ 

Since the imaginary part of the tadpole diagram vanishes, it does not contribute to 
the damping rate. For a fast moving particle with hard external momenta it is sufficient 
to use the effective propagator for one of the internal lines and bare propagators and 
vertices otherwise according to the resummation scheme of Braaten and Pisarski B , |l0| . 
Because of the Bose enhancement factor l+n# for the gluon loop, the diagram of Figjl](a) 
will give the leading order contribution to the imaginary part of the polarization tensor, 
whereas the quark and ghost loop can be neglected due to the Pauli blocking factor 1 — nj 
and will contribute only to higher order in the coupling constant. Therefore up to leading 
order we only have to calculate the three-gluon one- loop resummation graph of Fig|l|(a). 

In the static limit the contributions from the hard thermal loops for hot QCD are 

1,1,0] 

U L (k ^0,K)=n 00 (k ^0,K) = m 2 el = g 2 T 2 h + ^ + JL^^j (8) 



for the longitudinal part and 



n T (A:o^O,^) =0 (9) 



for the transverse part of the gluon polarization tensor. Substituting Eqs.(||),(|]) into 
Eq.(9), one gets 

Ai = At = h (10) 

The resummation method for perturbation theory at finite temperature as proposed 
by Braaten and Pisarski, which resums HTL contributions, works well if the momentum 
scale is gT. But in non-abelian gauge theories there is another pernicious barrier for per- 
turbation theory besides the weak coupling limit: Static magnetic fields are not screened 
at the HTL level and this may lead to a breakdown of perturbation theory at high loop 
order. The corresponding infrared singularities are usually expected to be cured by the 
dynamical generation of a magnetic screening mass m g (g 2 T) but the nature of 

the latter is still unclear. 

It is known that damping of collective modes in a QGP occurs. This dynamical 
screening, as opposed to static (Debye) screening, makes it reasonable to consider a 
damping rate for gluons, that is to assign them a finite width when studying the collective 
behaviour of a QGP. Assuming that the damping rate of the transverse gluon is z/p, we 



can write the transverse gluon propagator as 



so that the total resummed gluon propagator becomes 

G'^(k) = A^ u (k)M T (k) + B^(k)A L (k) + D^(k)^{k). (12) 
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The bare propagator and vertex are given by 



A^(k) 

(p, q, k) 



1 



9v> - (1 - 



kfj, k v 



k 2 yypu s/ fc2 

(P - q) P 9pu + (q~ k)^ g vp + (k - p) v g^ p . 



(13) 
(14) 



By making use of Eqs. (pl]) -(|l4T), one can evaluate the diagram of Fig.|](a) and obtain for 
F from Eq.(6) after a rather lengthy calculation 



9 2 N C 



d D q 4|P| 



(2tt)D {p + qf 



A' T {q)-A L {q) -{!-£) 



[P ■ Qf 



(15) 



where the hats on the three-vectors P and Q denote unit vectors. 

In the following we will apply the formalism of thermal field dynamics (TFD), since 
the imaginary part of the "1 — 2" component of the Greens function £12 is easier evaluated 
in TFD than in the imaginary time formalisme. By using the well known relation between 
the imaginary parts of retarded Greens functions and "1 — 2" components of Greens 
functions in TFD 0,0 

e P\po\ _ 1 

-ImSi2, 



ImS 



one immediately obtains 



ImF 



-9 2 N C 6 



P\po\ 



2e /3|po|/2 



d D q 
\D 



(16) 



2e /3|p |/2 J (2tt) 



A\P\ 2 A 12 {p + q) 



A 1 T 2 (q)-A 1 L 2 (q)-(l-0 (p.Q) 2 A 12 (q) 



The thermal propagators of 1 — 2 type have the form [14] 

n B (q)e-^ 2 , 



AP(<7) 
A?(q) 



2m5{q 2 - 
2 ImA T 



(17) 



(18) 



1 



1 



(q - iu T ) 2 ~ Q 2 



(q + iv T y - Q 2 



m 



9 j 



n B (q)e 



-P\qo\/2 



2m5{q 2 



■ n B (q)e 



-PM/2 



Ai 2 (?) 



2Tri8{q 2 ) ■n B {q)e 13 ^ 2 . 



(19) 
(20) 



We included the magnetic mass m g in the propagators, since both self-consistent damping 
and magnetic screening are non-perturbative effects of order g 2 T so that one cannot 
neglect a priori one effect compared to the other. In addition, including m g will enable 
us to compare our self-consistent results with the previous results, where the infrared 
divergence was regularized by a magnetic mass m g alone. 

Eq.(|l8|) contains the Bose-Einstein distribution function n B , which for soft momen- 
tum, q <C p ~ T, reduces to 



1 



oP\qo\ 



- — r, n B (p + q) ~n B {p). 



(21) 
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Substituting Eqs,(|l8| ) -(|2l| ) into Eq.(|17|), one can evaluate the integration by dimensional 
regularization at finite temperature p5|. The final result is 



ImF = - g2A : r|P| ( In f d 9 + 1 - - arctan ^ + 0(g 2 ) ) . (22) 
Air \ m 2 g + v?p 7r v T ' J 

From the definition Eq.(^) one obtains 

ImF g 2 N c T I ' ml, 2 m g \ 

VT = -^\=-^\ ln 2? 2 +1- -arctan^ . 23 

2|.P| 8-7T \ m~ + 7T I^T / 




We can rewrite Eq. ( |23D in terms of a new energy scale defined by m g with dimension- 
less variables vt = vx/nig, T = T/m g , \xt = /Ur/m 9 and 

as 

5 2 -/V c T / , . , .,. 2 1 \ 

z/y = 2 In m e z — m(l + u T ) + 1 arctan 3— . (25) 

87T V 7T vt J 

This is a self-consistent equation for the rescaled damping rate Vt and it is gauge inde- 
pendent to leading order in the coupling constant. 

If one neglects the finite width effects for the internal particles (z/j 1 1) Eq,(^5|) 
yields 

_ g 2 N c f _ 2 

»T= g7r lnm ei , (26) 



which recovers the previous results from Refs. |5[ , J/]] , |T( 

If on the other hand selfconsistent damping dominates magnetic screening (z/j- 2> 1) 
Eq.(p5p reduces to 

g 2 N c f / m 2 e X 



u T = ±-^[\ nT f + l , (27) 



87T \ Z/j, 



which can be rewritten as 



g 2 N c T f m 2 
8vr 1 n za 2 



^ T = ^- [ln-f + 11 , (28) 



where m g drops out completely. This again is a selfconsistent equation for ut and the 
infrared divergences encountered in previous calculations are removed self-consistently. 

The selfconsistent Eqs. (|25|) and (^) will prove to have very similar solutions for 
T > 1, but differ from the solutions for Eq. (|26|). 

One can calculate the damping rate for a fast quark going through a QGP at finite 
temperature and chemical potential in a way similar to the preceeding calculation for 
gluons by evaluating 

1 

Let us consider first the quark self-energy shown in Fig|2|, where we have included 
screening and damping rate effects by means of the effective gluon propagator ||. For 



r ImTr(^)| po=E . (29) 
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hard quarks {{p q ) ~ T) it is sufficient to use the bare quark propagator and bare vertices 
in the hard thermal loop approximation. 

We will calculate the imaginary part of the self-energy at finite temperature and 
chemical potential using TFD. For the hard particles one can show that the main contri- 
bution to the interaction rate comes from the soft momentum transfer range, i.e. q ~ gT 
in the weak coupling limit. 

Applying the usual Feynman rules of QCD to Fig.||, one obtains 

d D q Tvi>^{i> + i)Y 



SO, J 



(2*) 



D 



(p + qf 



-G 



/J.U1 



(30) 



with Ct 



N 2 -l 



jf = - 2N - the Casimir for SU(N). By using Eqs.(Q)-(|5|), Eq.pOj) can be cast into 
the form 



H 



-9 2 C 



d D q 2\P\ 2 
(2vr)^ (p + q) 2 



A' T (q)-A L (q)-(l-0 



(P ■ Qf 



(31) 



Then using the hot propagators in TFD for gluons as given by Eq.(^) and for quarks as 
follows 

1 



iA u (p) = -iA 22 (p) = i> 



2vr 



(0{po)e x v /2 n f {x p )+ 



jr — te 

e(-p )e-^ 2 n f (-x p ))8(p 2 ) 



iA\ 2 {p) = -iA 2 i(p) = 



e(p )e x ^ 2 n f (x p ) + 
e(-p )e-^ 2 n f (-x p ) 



5(p 2 



(32) 



where nf(x p ) = [e Xp + 1] 1 with x p = P(po + denotes the Fermi-Dirac distribution 
function, we obtain 



ImH 



-g 



2 ^/(0 e "^ 22|P|2sm2 ^( 2 -) 2 

A?(q) - A 12 (q) - (1 - Z)cos 2 9A 12 (q)) . 



(27T) 



1 



> ^5 [(p + q) 
1 



Q 2 -\- v Q 2 + m 2 



i-6 



cos 



Q 2 



(33) 



(34) 



Here we used the following relation between the retarded self-energy H and the 1-2 
component of the self-energy H\ 2 in TFD 

Ml 

ImH = — — lmH 12 , (35) 



where 



COS ( 



8111 ( 



sm2(f> po+fl 
0{po)e x v' 4 + 6\ 



-Po)e 



'po+P- 



\l q x v/^ -\- e x p/2 
9{po)e- x ^ A + 6{-p )e x v' A 

\J _|_ g—Xp/2 



and 



(36) 
(37) 
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Evaluating the integration in Eq.(^J) by dimensional regularization at finite temperature 
fi"5|| we get as the final result for the width of the quarks 



g 2 C f T„ mi, _ 2 m„ N 
v g = - — i— (In 5 el 9 + 1 - - arctan — £). 38 



Rewritting Eq.(38) as before in terms of dimensionless variables f g = v q /m g , v?, T, 
fh e i leads to 



9 2 Cff 
8vr 



2 1 \ 

2 In m e i — ln(l + p£) + 1 arctan — I , (39) 

7T i^t / 



which gives v g as a function of z/y and T, where vt is determined self-consistently through 

Eq. (Hi- 
Setting the finite width of the internal particle to zero, one recovers the result given 
in @. 

One thus finds that the damping rates for the quark and gluon have the same be- 
haviour with temperature and chemical potential to leading order in the coupling con- 
stant and differ only by a constant factor Cf/N c , which counts the number of processes 
available for color relaxation for quarks as compared to gluons. 



3 Numerical Analysis of Color Conductivity for a QGP 

In a QGP the interaction rates or frequencies among particles should be sensitive to 
variations of temperature and baryon density. Naively one might expect the interaction 
rates to increase with increasing temperature T or chemical potential /i, but one can see 
from the self-consistent Eqs.(p3|),(|3"8"|) that the interaction rates depend strongly only on 
temperature. The reason for this is that the chemical potential n. reflects only the baryon 
density in a QGP and does not describe the color charge density. 

We solved numerically the self-consistent Eqs.(p^),(|38|) for the gluon and quark damp- 
ing rates as functions of temperature for different values of the chemical potential in the 
weak coupling limit (a s = 0.1) for two quark flavors. 

Changing back from quantities u, T, where the energy scale is set by m g to ordinary 
temperature T and damping rate v we have to make an ansatz for the magnetic screening 
mass m g , which we know to be of order g 2 T: 

2 



m g = ca s T wither = — . (40) 



El 

The dependence of Eqs.(p3]),([3^) on c is given by a linear decrease of the asymptotic 
slope with increasing m g for fixed T, Aut / ' Ac\T= c <mst = —ttcx s T = const. For Eq.(p6|) on 
the other hand the change of the slope of vt as a function of c is given by Avt / Ac = — 3/c, 
so that for large c the asymptotic slope does practically not change. 

In the following we will adopt a value of c = 3 as estimated from lattice gauge theories 



1 17(1 and classical considerations fl8fl . 

For the considered values of T around 200 MeV in Fig. |3| we find an only slightly flater 
curve for Eq.(^3|), where the IR-singularities were regularized through a combination of 
selfconsistent damping and magnetic screening with m g = 3a s T compared to Eq.(p8|), 
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where selfconsistent damping alone was used. The latter case (limit !*r>l) corresponds 
to c = 0, that is vanishing m g in Eq.(H). 

For Eq.(p^) on the other hand, where the IR-singularities were regularized through 
a magnetic screening mass alone we find a steeper asymptotic slope. 

The results as given in Fig. || show that, as expected, the damping rate is not sensitive 
to the changes of the chemical potential, and only slightly enhanced with increasing [i. 
Otherwise it is approximately proportional to temperature. 

The influence of chemical potential is highest at very low temperature, corresponding 
to a cold (and dense) QGP, like for instance the core of a neutron star. For hot and dilute 
QGP as expected from ultrarelativistic heavy ion collisions the influence of chemical 
potential on the damping rate becomes negligible. 

From the transport theory of the QGP one can derive the static color conductivity 
as given in §],[§: 

u? 

^ = V' ^ = m e V3- ( 41 ) 

The color conductivity describes the ability of color transport in a QGP. It is defined 
as the ratio between the plasma frequency squared and the interaction rate T of the 
particles in a QGP, where T is the inverse color relaxation time and therefore given by 
the damping rate v. Since the interaction rate increases only linearly with temperature 
and logarithmically with chemical potential, but the plasma frequency squared increases 
quadratically with T as well as /U, color conductivity is enhanced with increasing T and 
H. 

As a numerical results for a one obtains the curves given in Fig.|| by using the self- 
consistent solutions for the damping rates. One finds that the quark chemical potential 
has some influence on the color conductivity, especially at low temperatures where our 
assumption of a finite chemical potential does apply. 



4 Discussion and Conclusions 

In summary, by taking into account the quark chemical potential and the transverse gluon 
damping effects, we obtain the self-consistent equation for the damping rate of gluons 
and quarks in real time finite temperature QCD language, where the infrared divergence 
is removed self-consistently and the magnetic screening mass leads to a redefinition of 
the energy scale. 

In addition we have solved the equation for the damping rate numerically and calcu- 
lated the color conductivity for a QGP for different chemical potentials with an ansatz for 
the magnetic screening mass m g . The results show that the chemical potential enhances 
the color conductivity somewhat but has practically no effect on the damping rate for 
gluons and quarks. 

Without considering a chemical potential or self-consistent damping we can recover 
the results previously derived by other groups. 

We conclude that the quark chemical potential has influence on color transport only at 
low temperatures and high densities. Except for low temperatures color conductivity rises 
linearly with T, where the slope of the curves depends on the mechanisme to regularize 
the IR divergences. Because of this T dependence dissipation can not be neglected in 
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hydrodynamic descriptions of the expansion phase of a hot QGP in ultrarelativistic heavy 
ion collisions. 
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5 Figure captions 



Figure 1: The resummed one-loop gluon self-energy graphs: (a) tadpole, (b) gluon loop, 
(c) fermion loop, (d) ghost loop. 



Figure 2: The resummed one-loop quark self-energy graph. 



Figure 3: Damping rate it for gluons versus temperature T for different quark chem- 
ical potentials, with a s = 0.1. Dashed lines: without selfconsistent damping for three 
different chemical potentials /j, = 0, 0.5, 1.0 GeV from bottom to top. Dotted lines: 
selfconsistent damping without magnetic screening, chemical potential \i as before. Solid 
lines: selfconsistent damping and magnetic screening, ji as before. 



Figure 4: Color conductivity a versus temperature T for different quark chemical poten- 
tials , with a s = 0.1. Coding of lines as in Fig. 3. 
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Fig. 1(b) 
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